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ABSTRACT 

We report time series and variability information for the optical identifications of X-ray sources in 47 Tucanae 
reported in Paper I (at least 22 CVs and 29 active binaries). The radial distribution of the CVs is indistinguishable 
from that of the MSPs detected by Freire et al. (2001). A study of the eight CVs with secure orbital periods (two 
obtained from the Chandra study of Grindlay et al. 2001a) shows that the 47 Tuc CVs have fainter accretion 
disks, in the V band, than field CVs with similar periods. These faint disks and the faint absolute magnitudes (My) 
of the 47 Tuc CVs suggests they have low accretion rates. One possible explanation is that the 47 Tuc objects 
may be a more representative sample of CVs, down to our detection threshold, than the CVs found in the field 
(where many low accretion rate systems are believed to be undiscovered), showing the advantages of deep globular 
cluster observations. The median Fx/Fopt value for the 47 Tuc CVs is higher than that of all known classes of 
field CV, partly because of the faint My values and partly because of the relatively high X-ray luminosities {L^). 
The latter are only seen in DQ Her systems in the field, but the 47 Tuc CVs are much fainter optically than most 
field DQ Hers. Previous work by Edmonds et al. (1999) has shown that the 4 brightest CVs in NGC 6397 have 
optical spectra and broadband colors that are consistent with DQ Hers having lower than average accretion rates. 
Some combination of magnetic behavior and low accretion rates may be able to explain our observations, but the 
results at present are ambiguous, since no class of field CV has distributions of both Lj and My that are consistent 
with those of the 47 Tuc CVs. 

The radial distribution of the X-ray detected active binaries is indistinguishable from that of the much larger 
sample of optical variables (eclipsing and contact binaries and BY Dra variables) detected in previous WFPC2 
studies by Albrow et al (2001). The X-ray properties of these objects (luminosity, hardness ratios and variability) 
are consistent with those of active binaries found in field studies, and the Fx/Fopt distribution is significantly 
different from those of the CVs and the millisecond pulsars (MSPs) that are detected (or possibly detected) in the 
optical. Despite these results, we examine the possibility that a few of the active binaries are MSPs with main 
sequence companions resulting from double exchanges in the crowded core of 47 Tuc. No solid evidence is found 
for a significant population of such objects, and therefore using the methods of Grindlay et al. (2002) we estimate 
that the number of MSPs in 47 Tuc with luminosities > lO-'^'erg s"' is ~30^0, near the previous lower limit. 

We present the results of a new, deeper search for faint low-mass X-ray binaries (LMXBs) in quiescence. One 
reasonable and one marginal candidate for optical identification of a quiescent LMXB was found (one is already 
known). Finally, it is shown that the periods of the blue variables showing little or no evidence for X-ray emission 
are too long for Roche lobe filling (if the variations are ellipsoidal). These blue variables also show no evidence 
for the large flickering levels seen in comparably bright CVs. At present we have no satisfactory explanation for 
these objects, although some may be detached WD-main sequence star binaries. 

Subject headings: binaries: general - globular clusters: individual (47 Tucanae) - techniques: photometric - 
X-rays: binaries 



L INTRODUCTION 

Binaries are well known to have a profound impact on the 
dynamical evolution of globular clusters (Hut et al. 1992 and 
Edmonds et al. 2003, hereafter Paper I), and they offer an op- 
portunity to study the results of stellar interactions. They also 
offer the crucial advantages of studying binary systems at the 
same, well-determined distance, age, metallicity and redden- 
ing. Accurate distances are particularly important for reduc- 
ing uncertainties in luminosities, and making reliable cluster- 
to-cluster comparisons and cluster-to-field comparisons. 

Despite these promising goals, observational progress in de- 
tecting globular cluster binaries was initially slow, mainly be- 
cause of spatial resolution and sensitivity limitations. However, 
the superb imaging capabilities of Chandra and HST has now 
resulted in the detection of large numbers of compact binaries 



in several globular clusters (see Grindlay et al. 2001a, here- 
after GHEOla and references in Paper I). In particular, using 
the Chandra data of GHEOla and the HST data of Gilliland et 
al. (2000; //5r program GO-8267), a large population of X-ray 
binaries with optical identifications has been found in 47 Tuc, 
with over 50 binaries reported in Paper I. There, the details of 
the astrometry and photometry for this sample were reported, 
with most of the systems being identified as either accreting 
white dwarf-main sequence star binaries (cataclysmic variables 
or CVs) or chromospherically active binaries. 

This paper reports HST optical time series and other detailed 
analysis for the sample of binaries reported in Paper I. Variabil- 
ity for most of the CVs is detected, confirming the CV identifi- 
cations given in Paper I based on absolute photometry alone. In 
particular flickering for the brighter optical IDs is found, plus 
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periods for several of the higher inclination systems and some 
long-term variability. The CV periods and the absolute mag- 
nitudes determined in Paper I are compared with predictions 
for Roche-lobe filling secondaries and are compared with CVs 
found in the field. Time series for the active binaries are also 
presented here, with an emphasis on the variables not discussed 
by Albrow et al. (2001; hereafter AGBOl). 

The radial distribution (radial offsets from cluster center) of 
the CVs and active binaries discovered in Paper I are compared 
with the radial distribution of MSPs from Freire et al. (2001) 
and with the total stellar population, to test whether these dis- 
tributions are consistent with the expectations of mass segrega- 
tion. Comparisons between the radial distributions of the X-ray 
detected active binaries and the larger sample of binaries dis- 
covered by AGBOl are also made. 

Other important diagnostics examined here are the X-ray lu- 
minosities and optical magnitudes. The flux ratio between the 
two was shown by Richman (1996) and Verbunt et al. (1997; 
hereafter VBR97) to be inversely proportional to the CV ac- 
cretion rate (for non-magnetic systems), and here we make de- 
tailed comparisons between the 47 Tuc CVs and field CVs to 
help constrain the CV accretion rates in 47 Tuc. The X-ray to 
optical flux ratio is also useful in searching for additional qui- 
escent low-mass X-ray binaries (qLMXBs) besides X5 and X7 
(GHEOla and Heinke et al. 2003). We also discuss the X-ray 
luminosity and absolute magnitude distributions of the 47 Tuc 
CVs and compare to field systems. 

The time series will be described in § 2, followed by an anal- 
ysis section (§ 3), including a study of the radial distribution of 
the sources and their X-ray to optical flux ratios. This will be 
followed by the discussion in § 4, and an interpretation of the 
results given in both Papers I and II. 

2. TIME SERIES 

2.1. Optical data 

Paper I discusses the photometric results for the optical IDs 
(including the CMDs) in detail. Exquisite time series were pro- 
duced in the V and / bands for the GO-8267 dataset, and are 
analyzed in this paper. Detailed simulations were carried out 
by AGBOl, giving the recovery rate of simulated variables as 
a function of period, V -magnitude and signal amplitude for the 
PCI and WF2 chips. For amplitudes of A(intensity)/intensity 
= 0.1, the recovery rate for PCI was ~100% from flie MSTO 
to y=23.5 for all input periods below 3.2 days, and for am- 
plitudes of 0.01 the recovery rate fell from ^70% at y=20 to 
<10% at y=21.5 (AGBOl). The WF2 recovery rates extend to 
slightly fainter magnitudes (AGBOl), but the generally higher 
crowding and saturation levels in the WF2 images counteract 
this gain, for variable searches. 

2.2. Cataclysmic variables 

2.2 A. GO-8267 data 

For the CV candidates found in flie GO-8267 FoV (Paper I) 
we have analyzed the high quality time series produced by 
the planet and binary searches of Gilliland et al. (2000) and 
AGBOl. Time series were produced using difference image 
analysis combined with aperture and PSF fitting photometry 
(see GilHland et al. 2000 and AGBOl for more details). The 
size of the region used in the photometry extractions was de- 
creased to as small as one pixel in the cases where plausible 
X-ray IDs were located near bright stars. Least-squares fits to 



sinusoids and the Lomb-Scargle periodogram were then applied 
to search for periodic variations in the time series of the CV 
candidates. Seven certain periodic variables were confirmed in 
fliis way: Wlopt, WSopt, W15opt, W21opt, W34opt, AKO 9 and 
the marginal ID for W71. We discuss each of these here in turn. 

The ID for Wl is either a 2.89 hr or a 5.78 hr period vari- 
able (eUipsoidal variations). Examination of the light curves 
folded at the 2.89 and 5.78 hr periods, and of the power spec- 
tra, do not show unambiguous departures from sinusoidal light 
curves (Fig. 1), although the V-band phase plot does show some 
evidence for asymmetry at a little less than the 3-a level, argu- 
ing for the longer period. We can use equation (2.102) from 
Warner (1995) to place limits on the absolute magnitude of the 
secondary: 

Mvi2)=l6J-n.llogPorbih), (1) 

where My(2) is the absolute magnitude of the secondary and 
Porh is the orbital period measured in hours (where 2 hr< Port ^ 
10 hr). This relationship is appropriate for field CVs of roughly 
solar metalUcity. At fixed mass, stars with lower metallicity will 
have brighter secondaries (Stehle, Kolb & Ritter 1997), and so 
equation 1 gives a lower limit to the luminosity of secondary 
stars in 47 Tuc (an upper limit to My). For Wlopt (My = 5.86) 
using a 5 .78 hr period gives a Umit of My = 8.2 from equation 1 , 
while using a 2.89 hr period gives My = 11.6. The longer pe- 
riod is therefore favored unless the accretion disk completely 
dominates in the V-band. This may be possible if Wl is a 
non-magnetic nova-like CV, explaining the relatively blue V-I 
color of Wlopt. However, the Fx/Fopt value for this system is 
much more consistent with a magnetic CV or a DN system (see 
Fig. 19), and a period of 2.89 hrs would be shorter than any of 
the UX UMa systems given in Fig. 10a, and all but one of the 
DQ Her systems. Modehng of the accretion disk and secondary 
in this system may help determine which period is correct. 

Note that both W8opt and W15opt are optically crowded and 
the absolute photometry in V and / is relatively poor (an / mag- 
nitude could not be self-consistently derived for W15opt). How- 
ever, since this project was optimised for time series accuracy, 
the quality of the time series is much higher than that of the ab- 
solute photometry. Both W8opt and W15opt are clearly eclipsing 
systems (see Fig. 2), with unambiguous period determinations 
of 2.86 hr (at the upper edge of the period gap for field CVs) 
and 4.23 hr respectively. Both of the associated X-ray sources 
for these CVs are very hard (GHEOla). Detection of these ob- 
jects as eclipsing strongly supports the argument of GHEOla 
that self-absorption by the accretion disks are responsible for 
the hard spectra. See § 2.2.3 for more details about these two 
binaries. 

The ID for W21 is very near to a strong saturation trail and 
evidence for a periodic signal was only detected in the single- 
pixel analysis. Fig. 3 shows the secure 52 minute variation (this 
may also be a 104 minute orbital period). The lack of strong 
signal at 104 minutes in the power spectrum and the large am- 
plitude of the time series (0.28±0.14) tentatively favor 52 min- 
utes being the orbital period. Poor phase coverage at phases 
0.4-0.75 makes it difficult to distinguish between these two pe- 
riods, although the phase plot does show suggestive evidence 
of asymmetry. In the 52 minute period case the secondary is 
possibly a degenerate star since 52 minutes is well below the 
canonical period minimum for field CVs of ~ 75 min, for non- 
degenerate secondaries. However, some field CVs with non- 
degenerate secondaries do have periods about this smaU (e.g. a 
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59 minute orbital period for V485 Cen; Augusteijn et al. 1996). 

The light curve for W34opt is discussed in detail in Edmonds 
et al. (2002b). Although the strong 97.45 minute period is very 
close to the HST orbital period (96.5 minutes), Edmonds et al. 
(2002b) argue that this is not the result of an artifact showing up 
as a false signal at the HST orbit. The binary W36opt (AKO 9) is 
a long period eclipsing CV with a subgiant secondary and a pe- 
riod unambiguously determined to be 1.108 days (see also Ed- 
monds et al. 1996, GHEOla and Knigge et al. 2002). The blue 
variable and possible counterpart to W71 has a 1.19 hr power 
spectrum signal with a likely 2.37 hr orbital period, based on 
the small but significant asymmetry in the phase plot using the 
longer period (see Fig. 4). If this is a CV it would fall within 
the period gap for field CVs (2.2 < P^rb < 2.8hr). Although 
this variable may be a magnetic CV, and the source is relatively 
soft, consistent with an AM Her system, we favor the explana- 
tion (supported by the astrometry), that the Chandra source is 
not associated with the nearby blue, variable star (see § 3.4 for 
possible explanations of its nature). 

For the other CV candidates, deeper searches for periodic 
variations were made and several marginal signals were found 
(see Fig. 5). For VI, the strongest V-band power spectrum peak 
in a reasonable period range corresponds to an orbital period of 
~7 hr (for elUpsoidal variations). For V2 peaks are found at 
13.67/2 hr (V) and 6.01/2 hr (/) and a strong signal is also seen 
in V at 1.44 days and a 1.5 day signal in I. The shortest of 
these corresponds to a plausible orbital period, while the two 
longest ones are unlikely to be orbital periods for this apparent 
MS secondary, unless it is evolved. For V3, peaks are found in 
both V and / corresponding to a possible orbital period of ~3.7 
hr (close to the 3.83 hr period present in the X-ray light curve; 
GHEOla). The ID for W45 shows a significant peak at 1.87 
hr in V (possibly a 3.75 hr orbital period). A weaker signal is 
found in / at the same period and a strong signal at 1.78 hr. 

The ID for W2 shows possible periods at 2.2, 5.9 and 8.2 hr. 
These differ from the statistically significant X-ray period of 
6.287 hr (see Fig. 6 for an X-ray phase plot and Fig. 7a for a V- 
band power spectrum and phase plot). The power spectrum is 
extremely noisy, at least partly because the counterpart lies on a 
diffraction spike from a nearby giant star. Flickering may be an 
extra source of noise. The ID for W120 shows a marginally 
significant peak in the V power spectrum of 1.32 hr, with a 
false-alarm probability (FAP) of 1.4x10"^ (see Fig. 7b). The 
power spectrum is otherwise very clean, with no other peaks 
having a FAP < 0.9. The possible counterpart to WIO has a 2.2 
or 4.4 hr period and may also be a CV, but crowding prevents 
us obtaining useful Umits on the color. 

Figure 8 shows plots of time series rms vs V and / for the 
PCI, WF2 and WF4 chips (the few WF3 variables are shown 
on WF4). In the cases where / values are not available (W15opt, 
W33opt, W45opt and W70opt) we have assumed that the object 
is found on the MS in the V vs V- / CMD. As expected, sig- 
nificantly higher rms values than average are found for stars 
with large amplitude, periodic variations (e.g. Wlopt, W8opt and 
W15opt). Among the CV candidates, significant non-periodic 
variations were found for VI, V2, V3 and W25opt. These ob- 
jects are clearly visible in the V and / images, and are not af- 
fected by diffraction spikes or saturation trails, so the source 
of the enhanced noise for these objects is probably flickering 
from the hot spot or accretion disk. Strong temporal changes 
in the PSFs correlated with telescope breathing caused likely 
spurious signals to appear in the power spectra near the HST 



orbital frequency (0.1725 mHz) for these four large amplitude 
variables. 

The same explanation cannot be confidently appUed to the 
possible enhanced noise for W33opt, W45opt, W70opt and 
W120opt. These objects are all so crowded that none of them 
were detected in the original V and / band images by AGBOl, 
and the enhanced noise could have some contribution from this 
crowding. Strong signals near the HST period are only found 
for W33opt in the / band. By contrast, the ID for W44opt is rel- 
atively isolated, though faint, and the higher rms values found 
here are also suggestive of non-periodic variations. 

2.2.2. GO-7503data 

Significant variability was found for several of the CV candi- 
dates by comparing the 3rd epoch images of the HST program 
GO-7503 with those obtained in the first two epochs (epoch 1: 
GO-5912, and epoch 2: GO-6467). Fig. 9 shows those cases 
where F300W variability was discovered (small amplitude vari- 
ability was also found for V2, a known DN, but is not shown 
here). Most of the variations are obvious from Fig. 9 with the 
exception of W44opt which is 0.36 mag fainter in the 1st epoch 
than it is in the 2nd (and 0.30 mag fainter in the 1st than in the 
3rd), and W122opt which is 0.87 mag fainter in the 1st epoch 
than the 2nd and 0.52 mag brighter in the 2nd epoch than the 
3rd. Some of these variations may represent dwarf nova (DN) 
outbursts, especially W5 lopt and W56opt. A detailed analysis of 
the long-term variability of the CVs will be deferred to a later 
publication (where these results will be combined with subse- 
quent i/5r observations made with ACS). 

2.2.3. Period-Mv plot 

Having a sample of CVs and CV candidates with period and 
My determinations at a well-known distance allows us to exam- 
ine the period luminosity relationship expected for CVs. Figure 
10a shows a plot of period vs My for CVs and CV candidates in 
47 Tuc along with those of CVl and CV6 in NGC 6397 (Grind- 
lay et al. 2001b, hereafter GHEOlb; Kaluzny & Thompson 
2002), the DN in M5 (Neill et al. 2002), two CVs in NGC 6752 
(Bailyn et al. 1996) and the magnetic CV in M67 (GilUland et 
al. 1991). The straight line is equation 1 from Warner (1995), 
for solar-metallicity stars. This will shift to brighter magnitudes 
(by about 1 mag) for metal-poor secondaries, so CVs that are 
near or above this line are likely to have evolved secondaries. 
Some of the V-band fight will be from the accretion disk and 
hot spot (as well as the secondary) and this will also push the 
V magnitudes to brighter values (to the left in 10a). Figure 10b 
shows the plot for the blue variables described earlier. 

The Ukely CV W2opt and flie CV candidate V3 both lie close 
to or above the MS ridge-line in the V vs V-/ CMD (see the 
CMD for PCI in Paper I), showing that the secondary domi- 
nates the optical light. As expected these objects lie close to 
the Warner relationship. The secondary for W2opt might be 
evolved, but better photometry is needed to make this conclu- 
sion secure. The CVs VI and Wlopt fall below equation 1 be- 
cause the accretion disk and hot spot (and possibly the WD) 
make significant contributions to the light (the CMDs in Paper 1 
show that the y- / colors of Wlopt and VI are well to the blue 
of the MS). The light curve of the eclipsing CV W15opt (Fig. 2) 
implies that the accretion disk and hot spot are responsible for 
~ 60% of the V-band light, and increasing V by ^1 mag to iso- 
late the MS component brings W15opt into excellent agreement 
with equation 1. For W8opt around 40% of the V-band light ap- 
pears to come from a source other than the secondary. Given 
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the similarity of the hght contributions from the secondaries, 
equation 1 implies that the W8opt secondary should be about 2 
mags fainter than the WlSopt secondary, but instead it is almost 
1 mag brighter. Possible explanations are that: (1) we are see- 
ing a grazing eclipse with the inclination slightly less than 90 
degrees, (2) a faint MS star is contaminating the hght from the 
CV, or (3) the normalization for W8opt is incorrect because of 
the crowding. Examination of Fig. 2 shows that the elUpsoidal 
variations in WlSopt are close to the maximum value expected 
(see Russell 1945), but are around a factor of 2 smaller in W8opt, 
consistent with any of these three explanations. Similar (but 
less extreme) issues exist for the two NGC 6752 CVs (both of 
these systems He close to the MS in the V vs V- / CMD). The 
good agreement of the M5 DN with equation 1 implies that the 
V-band light must be dominated by the secondary. 

The variable W34opt has colors that are consistent with those 
of a CV {U-V = -0.12; y-/=1.03; Mv=9.3) and appropriately 
falls below equation 1. However, as noted by Edmonds et al. 
(2002b), the relatively large amplitude variation of this object 
in both V and / makes an MSP another possibility, since the 
light curve is similar to that of W29opt(47 Tuc W). Either possi- 
bility is consistent with the X-ray properties of this object. The 
variable W21opt (U-V=0.15; V -1=0.21; Mv=7.44) is also dif- 
ficult to categorize. Assuming that the orbital period is twice 
the observed period, the star still falls well below equation 1. 
Could an accretion disk for a CV with a period of 1.73 hr and 
an expected My (2)= 14 be this bright? Figure 9.8 of Warner 
(1995) shows that the 13 DNe with periods <2 hr mostly have 
My(disk)= 8-9, with 3 having Myidisk) < 7.1. However, we 
would expect a very blue U-V color if the disk is this bright. 
Alternatively, the binary could be a WD-WD CV or a WD-NS 
system in quiescence. Finally, AKO 9 with its subgiant sec- 
ondary is well outside the range of applicability of equation 1 . 

Figure 1 1 shows a period vs My plot (similar to Fig. 10) for 
the field CVs with estimated distances in the ROSAT study by 
VBR97. The orbital periods have been taken from Ritter & 
Kolb (1998). The nova-like CVs (systems that are magnetic 
or do not show outbursts) and the DNe (outbursting systems) 
have been separated into two plots to show systematic differ- 
ences between these two general classes of CV. The labeling of 
VBR97 for the different CV types is used. The nova-like CVs 
are generally further away from equation 1 than the DNe, show- 
ing that the relative hght contributions from the accretion disk 
and hot spot (compared to the secondary) are larger than for the 
DNe. The 47 Tuc CVs (and the other cluster CVs) appear to be 
more consistent with the DNe from VBR97 than the nova-like 
systems. 

By calculating the magnitude predicted by equation 1 (us- 
ing the binary period) and subtracting My, we quantify the off- 
sets from equation 1 and examine the light contribution from 
sources other than the secondary (mainly the accretion disk, hot 
spot and WD). This procedure will overestimate the contribu- 
tion from the disk and other non-secondary components, since 
equation 1 underestimates the contribution from a metal-poor 
secondary. Figure 12 shows the cumulative distributions of the 
offsets from equation 1 for the 47 Tuc and the field CVs. We 
have eUminated binaries with periods >10 hr to remain con- 
sistent with the period limitations on equation 1. Figure 12a 
shows all such systems and Fig. 12b shows only systems with 
periods above the CV period gap (the latter figure avoids the 
shortest period systems likely to have a greater light contamina- 
tion from the WD). For systems above the period gap, the offset 



distributions for the 47 Tuc CVs and the nova-like systems are 
distinguishable at the 99.55% level using the KS-test. A subset 
of the nova-like systems, the DQ Hers, are distinguishable from 
the 47 Tuc CVs at the 93.8% level (the nova-like and DQ Her 
offsets are distinguishable from each other at the 1.2% level). 
By contrast, the 47 Tuc CV offsets are distinguishable from the 
DNe offsets and the U Gem systems at only the 1.4% and 0.9% 
levels respectively. This suggests that the 47 Tuc CVs have 
accretion luminosities that are even lower than those found in 
U Gem systems and DNe since, as mentioned above, the non- 
secondary component in the 47 Tuc CVs is being underesti- 
mated here. This then suggests that the 47 Tuc CVs have lower 
accretion rates. More evidence supporting this conclusion is 
given in § 3.2, and discussion is given in § 4.1. 

2.3. Blue Variables 

As described in Paper I, several blue, variable stars discov- 
ered by AGBOl are not obviously associated with X-ray sources 
and have colors that are different from most cluster and field 
CVs. The time series for these blue variables are shown in 
Fig. 13 and the period versus My plot is shown in Fig. 10b. 
The variable 1V36 has a sinusoidal, 9.55 hr variation; 2V08 
and 2V30 have 4.8 hr and 28.3 hr variations, and 3V06, 3V07 
and 4V05 have 2.78 hr, 3.79 hr and 3.11 hr variations respec- 
tively. Note that the absolute magnitudes of the secondaries (if 
MS members) are likely to be significantly fainter than the ob- 
servedMv because these stars are all blue in V vs V-I. We have 
assumed that the orbital period is twice the detected period, as 
appropriate for most CVs (Roche-lobe filling, MS secondaries 
with ellipsoidal variations) with periodic variations detected in 
the optical. If this assumption is correct, then 1V36, 2V08, 
2 V30 and 3 V06 are unlikely to be filling their Roche lobes since 
they fall above equation 1 in Fig. 10b, and are therefore unhkely 
to be CVs (as akeady noted for 1 V36 by AGBOl). If any of the 
secondaries are evolved then they could fill their Roche lobe at 
the observed period, but such secondaries should be unusually 
red, and therefore to explain the observed blue V-I colors we 
probably require an unusually UV-bright, luminous component. 
Even if evolved companions were expected to give unusually 
bright disks, only 2V08 and 2V30 have U-V colors as blue 
as most of the other CVs, so the blue variables do not seem to 
have unusually bright disks. The blue color for 3V07 implies 
that only a small fraction of the optical light comes from a MS 
star, and therefore the secondary for this system is also Ukely to 
fall above the Warner relationship (the My offset of 3V07 from 
equation 1 in Fig. 10 of ~1 mag roughly equals its offset from 
the MS in the V\sV-I CMD. 

Another possibiUty is that the orbital periods are the same as 
the observed periods, and that the variations are caused by heat- 
ing of one side of an MS star In this case some of the blue vari- 
ables (like 1 V36) may not have to contain evolved secondaries 
for Roche lobe filling to be taking place. There are at least two 
problems with this explanation. Firstly, a UV-bright, relatively 
luminous component is probably required to cause measurable 
heating effects, but only 2V08 and 2V30 have U-V colors that 
are comparable to those of the bluest CVs. Secondly, because 
sinusoidal signals with low amplitudes (either from heating or 
ellipsoidal variations) are generally not detectable even in the 
brightest CVs because of noise from flickering. If the blue 
variables are CVs then we would expect flickering to mask the 
small sinusoidal signals that are observed. Instead, the rms val- 
ues for the blue variables are quite low (Fig. 8) implying low 
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amplitude flickering, despite having bluer colors and (presum- 
ably) relatively bright disks, if they are CVs. 

We note that these results are unlikely to be some sort of pho- 
tometric artifact (such as the effects of bad pixels), since several 
of these objccls exist and none of them show unusual PSF fits, 
or excessive levels of crowding. Also, the power spectra show 
little if any extra power near the HST orbital frequency. See 
Sect. 3.4 for more discussion about these imusual objects. 

2.4. Active binaries 

The phase plots for W12opt, W92opt, W137opt and W182opt 
(eclipsing binaries) and W41opt, W47opt and W163opt (contact 
binaries) are shown in AGBOl. Light curves for most of the 
BY Dras and red stragglers from AGBOl that have also been 
identified as active binary candidates are shown in Fig. 14, and 
phase plots for a selection of active binary candidates not found 
by AGBOl are shown in Fig. 15. Note that ambiguity exists in 
the period determination for several of the variables shown in 
Fig. 15 (especially for W22opt and W94opt). In these cases we 
have assumed that the orbital period is twice the observed pe- 
riod. Clearly, W38opt and W94opt are eclipsing binaries, and the 
other new objects are a likely mixture of BY Dra and W UMa 
systems (W66opt was classified by AGBOl as a non-ecUpsing 
contact binary). 

3. ANALYSIS 

3.1. Radial Distribution of Sources: Comparison Between 
Source Classes 

Here, we compare the radial offsets, with respect to the center 
of 47 Tuc ('radial distributions') of the different classes of bi- 
nary. Although the GO-7503 dataset has uncovered a number of 
likely CVs and at least 2 active binary candidates, for radial dis- 
tribution analysis we restrict ourselves to the GO-8267 dataset 
because it contains a much more complete sample of active bi- 
nary candidates. Figure 16 shows the radial distributions of: (1) 
the radio detected MSPs (plus W34), (2) the CVs, (3) the bina- 
ries discovered by AGBOl (with the relatively small fraction 
of X-ray bright systems and the two confirmed CVs removed), 
(4) the candidate active binaries or 'xABs' (MS or red straggler 
binaries either from AGBOl or this work with apparent X-ray 
counterparts) and (5) the blue variables. Note the similarity be- 
tween the CV and the MSP radial distributions (KS-test prob- 
abihty = 14.0%). This similarity is not surprising because the 
single MSPs should have masses of ~ 1 .4M0 and the MSP bi- 
naries only ^O.O2-O.2M0 more. The detected CV candidates 
should have masses ranging from '^O.95M0 for faint objects 
like V3 and W120opt (--0.4Mq secondary and ~O.55M0 pri- 
mary) up to ~1.6 M0 for AKO 9 (^O.8M0primary and sec- 
ondary). The 7 blue variables are sHghtly less centrally con- 
centrated than the CVs (KS-test probability = 88.2%) and the 
MSPs (KS-test probability = 93.8%). EHminating 2V08 and 
2V30, the two variables with colors most like white dwarfs or 
CVs, causes the above KS-test probabiUties to drop to 53% and 
55%, suggesting that the blue variables are heavy objects (2V08 
and 2V30 are relatively distant). However, the remaining sam- 
ple of 5 is very small for statistical tests. 

The radial distribution of the AGBOl binaries is slightly less 
centrally concentrated than those of the CVs and MSPs (KS-test 
probabilities of 77% and 81% respectively). These results are 
only marginally statistically significant, but are a likely result 
of lower masses for some of the AGBOl binaries, especially 
the faintest ones. The active binaries appear to be marginally 



more centrally concentrated than the AGBOl binaries, but the 
KS-test probability (45%) is not statistically significant. If real, 
this difference likely results from the generally brighter V mag- 
nitudes (and higher masses) of the active binaries compared to 
the AGBOl binaries (see Fig 16b), since the V distributions dif- 
fer at the 90% probability level. Also, the V magnitudes of the 
X-ray bright BY Dras are brighter than the V magnitudes of 
the much larger BY Dra population from AGBOl (see Fig 16d; 
KS test probabihty = 93.5%). The active binary periods ex- 
tend to higher values than the periods of the CVs and the MSPs 
(Fig. 16c), and most strikingly they are significantly shorter 
than those of the AGBOl binaries (KS probabihty = 98.9%). 

We now consider the AGBOl binaries and active binaries 
with y < 19 and V > 19 (this divides the AGBOl binaries into 
two evenly sized groups). In the bright subsample the radial 
distributions of the AGBOl binaries and the active binaries are 
indistinguishable from each other (KS-test probability = 16%). 
The radial distributions of both the AGBOl and the bright active 
binaries are indistinguishable from those of the MSPs (KS-test 
probabihties = 30% and 29%) and CVs (KS-test probabiUties = 
16% and 9%; see Fig. 17a). 

Considering the fainter binaries (V > 19), the radial distribu- 
tion of the faint AGBOl binaries and the faint active binaries are 
indistinguishable from each other (KS-test probabihty = 12%). 
The radial distributions of both the faint AGBOl binaries and 
the faint active binaries (see Figure 17b) are less centrally con- 
centrated than those of the MSPs (KS-test probabilities = 97.0% 
and 83%), with similar results for the CVs (KS-test probabili- 
ties = 93% and 85%). These results are consistent with the 
faint active binaries being less massive than the bright active 
binaries. 

Figure 17c shows the radial distribution of the sources in the 
GO-8267 FoV that have not been identified, plotted with the 
bright and faint active binaries and the general stellar popula- 
tion in three different V-band ranges. The unidentified sources 
are indistinguishable from the bright active binaries, the CVs 
and the MSPs, showing that they are consistent with any of 
these groups, or a combination of all three types. The faint 
active binaries have a radial distribution that is almost identical 
to that of the stars with 17 < V < 18 (probability = 0.73%), and 
is considerably more concentrated than the stars with 19 < V < 
22 (probability = 84.7%). All of these results are consistent 
with expectations for mass segregation of binaries and stars in 
47 Tuc, and with the statement of AGBOl that incompleteness 
as a function of radius and magnitude is expected to be small 
even for stars as faint as V^21.5. 

The most striking result is that the period distribution of the 
faint active binaries differs from that of the faint AGBOl bi- 
naries at the 99.99% probability level. Figure 17d shows that 
while all of the faint active binaries have periods <0.7 days, 
only ~30% of the AGBOl binaries have periods below this 
value. We discuss in § 4.2 how these results imply the pres- 
ence of a stronger relationship between luminosity and period 
than found in field objects. 

3.2. X-ray to optical flux ratio 

3.2.1. Comparison Between Source Classes 

A useful, distance-independent diagnostic for stars detected 
in both the optical and X-ray wavelength bands is the ratio of 
the X-ray flux to the optical flux. Figures 18a and 18b show 
plots of this ratio vs My for nova-like CVs and DNe in the 
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field (VBR97). The values in the 0.5-2.5 keV band from 
VBR97 were scaled to be consistent with the shghtly absorbed, 
1 keV thermal bremmstrahlung spectra (0.5-2.5 keV) assumed 
by GHEOla, and used here. The Fopt values were derived from 
Fopt = lO"''-'*^"^ "^^. In cases where the distances to the field CVs 
of VBR97 were unknown, the systems were plotted on the left 
side of Figures 18a and 18b. 

The majority of field DNe have Fx/Fgpi > 0.01 with the 
smallest ratios being dominated by Z Cam systems (with rela- 
tively high accretion rates). The remaining DN classes (U Gem 
and SU UMa systems) mostly have Fx I Fopt > 0.1. Among the 
nova-like CVs the UX UMa systems have FxjFgpt values or 
limits mostly around 0.01. These ranges show a clear trend that 
higher FxlF„p, values correspond to lower accretion rates for 
non-magnetic CVs, as pointed out by VBR97 for this sample of 
field CVs, and by Richman (1996). Most magnetic CVs have 
Fx I Fopt > 0. 1 , although it is unclear if Fx /Fopt is anti-correlated 
with accretion rate for these systems. 

Fig. 19a shows the FxjFgpt values for the 47 Tuc CV can- 
didates and MSPs (plus other blue IDs) and Fig. 19b shows 
the Fx/Fopi values for the active binary candidates. Limits are 
shown where appropriate. Note the clear differences between 
Figures 19a and 19b, with generally much larger Fx /Fopt val- 
ues for the CVs and the MSPs. With the exception of AKO 9, 
the Fx /Fopt values for 47 Tuc CVs and CV candidates are all 
> 0.1, suggesting that the population of 47 Tuc CVs, if dom- 
inated by non-magnetic systems, are low accretion rate CVs 
(U Gem and SU UMa systems, where the latter sub-class are 
beUeved to have the lowest accretion rates among field CVs). 
Besides V3, the likely DQ Her system VI (GHEOla) has the 
highest Fx /Fop, value. 

Cumulative distributions of Fx / Fopt for the 47 Tuc CVs, ac- 
tive binaries and MSPs are compared in Fig. 20a. Only the 
objects detected in V have been included rather than those with 
limits, with the exception of 47 Tuc T. The clear differences be- 
tween the CV and the active binary distributions (and between 
the MSP and active binary distributions) are confirmed by the 
median values given in Table 1 and the large KS-test probabil- 
ities shown in Table 2. The MSP distribution only contains 4 
systems and the other binary MSPs detected in the radio fall 
below our optical detection threshold (but are mostly detected 
by Chandra), so the true MSP distribution will extend to higher 
Fx /Fopt values. The only region of significant overlap between 
the active binaries and either the CVs or MSPs is found with 
0.1 < Fx /Fopt < 0.5. All of these apparently extreme active 
binaries are found in the faint group and have radial distribu- 
tions that are less centrally concentrated than either the CVs, 
the MSPs, or the bright active binaries. 

3.2.2. Cataclysmic variables 

In Fig. 20b the 47 Tuc Fx /Fopt distribution for the CVs is 
compared with those found for field CVs using the data given in 
VBR97. Table 1 shows the median values for each distribution 
and Table 2 compares the KS-test probabilities. The trend of in- 
creasing Fx/Fopi for decreasing accretion rate (Richman 1996 
& VBR97) is clear, with the UX UMa (nova-hke) systems and 
the SU UMa (DNe) systems at opposite ends of the range, and 
the Z Cam and U Gem systems falling in between. The 47 Tuc 
Fx /Fopt values are higher on average than all of the various CV 
subgroups, including 3 old novae and 3 double-degenerate sys- 
tems not shown in Tables 1 and 2. In a separate calculation we 
have also included the field systems with X-ray limits (rather 



than just detections) given by VBR97 and this only increases 
the differences between the Fx/Fgpt values of the 47 Tuc CVs 
and the field CVs. 

To further investigate the origin of differences in Fx / Fopt be- 
tween the 47 Tuc CVs and the field CVs, we plot in Fig. 21a the 
cumulative distributions of the 47 Tuc CVs and field CVs and 
a similar plot in Fig. 21bforMy. Limits in L;^ (corresponding to 
X-ray detection Umits, not distance limits) have been included 
to increase the sample sizes of the field CVs, since only a subset 
of the CVs have estimated distances. The luminosities for most 
of the field CV classes are indistinguishable from each other, 
with only the DQ Her systems having a significantly different 
(more luminous) distribution, as pointed out by VBR97. The 
Lj distribution of the 47 Tuc CVs is distinguishable from the 
AM Her, DQ Her, UX UMa, Z Cam, U Gem and SU UMa sys- 
tems at the 97.6%, 22.4%, 97.3%, 72.2%, 99.3% and 99.8% 
levels respectively, so that only the DQ Her systems are ob- 
viously a good match to the 47 Tuc CVs. This would argue 
against the identification of the 47 Tuc CVs with low accretion 
rate CVs like DNe, if the accretion rate is proportional to 
(for example Warner (1995) suggests Lace — SOL;^). However, 
see the discussion of the My distribution below. 

Note that the majority of field AM Her and DQ Her systems 
from VBR97 were X-ray selected (as were the 47 Tuc CVs), 
unlike other classes of CV, implying that the relatively high 
X-ray luminosities could be partially selection effects, perhaps 
weakening the classification of the 47 Tuc CVs with magnetic 
systems. For example, crowding in X-rays is likely to have 
prevented the detection of some systems with low X-ray lumi- 
nosities, and the true cluster distribution will therefore extend 
to lower luminosities. Robust limits on incompleteness will be 
presented in future publications. Incompleteness will also ap- 
ply with the field samples, especially the magnetic systems. For 
example, although only 2 out of 10 of the DQ Her systems in 
the field sample analysed here involve upper Umits, a signifi- 
cant number of faint DQ Hers (and other CV classes) may have 
been missed in field surveys. The four most luminous DQ Her 
systems are at relatively large distances (>400 pc), suggesting 
that such bright systems are relatively rare, as pointed out by 
VBR97. However, despite this incompleteness at low luminosi- 
ties, the similarity between the 47 Tuc and DQ Her distributions 
at high luminosities (above ~ 7-8 x 10^' erg s"') is noteworthy, 
and VI, the brightest CV in 47 Tuc, has already been shown to 
be a likely DQ Her system. 

The My distribution of the 47 Tuc CVs is distinguishable 
from the AM Her, DQ Her, UX UMa, Z Cam, U Gem and 
SU UMa systems at the 8.6%, 99.97%, 99.99%, 98.80%, 1.03% 
and 96.2% levels respectively. Here, only the AM Her and the 
U Gem systems are good matches for the 47 Tuc CVs. A com- 
parison of Figures 18 and 19 shows that our optical detection 
limit of V 23 (because of crowding; see § 4.5) is a significant 
source of incompleteness for faint DNe searches, and may pre- 
vent the detection of a significant fraction of SU UMa systems. 
Clearly, the 47 Tuc CVs are significantly fainter in the opti- 
cal than the nova-like CVs detected in the field, including the 
DQ Her systems. Since metal-poor secondaries are expected 
to be brighter than solar metalUcity secondaries, at fixed mass, 
the difference in accretion luminosity between the 47 Tuc CVs 
and field DQ Hers will be even more dramatic than suggested 
by Fig. 21b. Finally, distances are available for only 2 old no- 
vae and none are available for double degenerate systems (both 
of these classes of CV may be extremely faint and beyond our 
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detection limits). 

The faint My values for the 47 Tuc CVs suggest that they 
have low accretion rates, consistent with the penod/My analy- 
sis, but conflicting with the interpretation based on their X-ray 
luminosities. In fact, an examination of Fig. 21 shows that the 
47 Tuc CVs have combined and My distributions that are 
not consistent with any known class of field CV (see § 4.1 for 
discussion of this issue). 



3.3. qLMXB search 

As expected, the optical ID for X5AV58 (Edmonds et al. 
2002a), a qLMXB, has a relatively high Fx/Fopt value of 47.8 
(Fig. 19; see also Heinke et al. 2003). Interestingly, the CV 
candidate V3 has a very similar Fx/Fgp, value. This ratio is 
higher than that of all known field and 47 Tuc CVs, represent- 
ing evidence for V3 also being a qLMXB. The possibility that 
V3 is a binary containing a neutron star has already been dis- 
cussed by GHEOla who argue that V3 may be an enshrouded 
MSP. For comparison, source B, the possible qLMXB in NGC 
6652 has Fx/Fgp, -^9.0 when the Fgp, value given by Heinke, 
Edmonds, & Grindlay (2001) is appropriately converted, while 
limits for X7 and the qLMXB in NGC 6397 are > 166 (Heinke 
et al. 2003) and >64 respectively. However, we caution that 
the Chandra and GO-8267 data were not obtained simultane- 
ously, and since V3 is highly variable both optically (see Pa- 
per 1) and in X-rays (Verbunt & Hasinger 1998; GHEOla) the 
high Fx/F„p, value derived here could be explained by a low 
optical state during the GO-8267 observations, or a high state 
during the Chandra observations. 

The two bright {L^ ^ lO^^erg s"') sources and likely 
qLMXBs in 47 Tuc (X5 and X7) and the likely qLMXBs in 
NGC 6397 and NGC 6440 are all soft sources, using the X-ray 
color definition of Grindlay et al. (2002), namely xcolor=2.5 
log([0.5-1.5keV]/[1.5-6keV]). Therefore, we have searched 
for evidence that any of the fainter, soft X-ray sources in 47 Tuc 
may be qLMXBs (besides V3). The source AKO 9 is a very soft 
source but is more hkely a CV than a qLMXB because it has 
been observed to have at least 2 outbursts without the detec- 
tion by ROSAT or RXTE of a 10^^ erg s"^ source in 47 Tuc 
(associated with an LMXB in outburst). The only other rea- 
sonable candidate is W17, a 103 ct source that is a factor of 
3 fainter than the qLMXB in NGC 6397. This source is out- 
side the GO-7503 FoV and unfortunately lies in the middle of a 
large saturation trail in the GO-8267 V and / images. However, 
useful limits can be set from the deep and relatively clean U 
image. The nearest detectable stars to the X-ray position are at 
2.58, 2.91 and 3.04 ct, and these have U = 22.4, 18.3, & 21.1. 
The fainter two of these are comparable to the CVs W44opt and 
W45opt and the brighter one to a system intermediate between 
AKO 9 and VI. We set an upper limit to the detection of a 
star at the nominal W17 X-ray position (ie assuming zero error 
in transforming to HST coordinates) of [/ ~ 24. This is sug- 
gestive of a faint optical counterpart like X5 (Edmonds et al. 
2002a) or a faint upper limit as with the qLMXB in NGC 6397 
(GHEOlb). Other possibiUties, besides a qLMXB, are that the 
source is a soft CV (e.g. an AM Her system) or an active bi- 
nary. However, for the latter possibility this source would be a 
factor of ^3 brighter than any of the soft (and non-flaring) ac- 
tive binaries. Followup observations with ACS should be useful 
in searching for optical counterparts and determining the nature 
of this source. 



3.4. Blue Variables 

As noted in Paper I, 1V36 may be a faint X-ray source (L^ 
< lO'^'erg s"'), arguing in favor of this object being a CV (al- 
though Edmonds et al. 2001 and Edmonds et al 2002b both 
give examples of faint X-ray sources with blue, variable optical 
counterparts that are MSPs, not CVs). Much deeper Chandra 
observations are required (and have been taken) to confirm this 
marginal detection of an X-ray source, but there are a number 
of arguments against the CV identification: (1) the period of 
19.1 hr (assuming the observed variations are elhpsoidal) is too 
long (see Fig. 10). Ellipsoidal variations should be seen in this 
low noise time series if the system is a CV and therefore Roche- 
lobe filling (unless the inclination is very low); (2) the position 
of this star in the color-color plot is different from most known 
field and cluster CVs; (3) the upper limit on Fx/Fopt for 1V36 
is lower than that of all 47 Tuc CVs (see Fig. 19); (4) there 
is no suggestion of the large flickering present in other opti- 
cally bright CVs such as VI, V2 and W25opt. Note that (2) and 
(4) likely rule out the possibility that several UX UMa systems 
with very small Fx jFopt values might have slipped beneath our 
X-ray detection limit. We therefore believe that the compre- 
hensive data available for 1V36 does not support the claim by 
Knigge et al. (2002) tiiat 1V36 is a strong CV candidate and 
look forward to seeing what STIS/FUV spectroscopy and vari- 
ability says about this system. 

Other explanations for 1V36 fare little better. The variable 
appears to be too faint and blue in Y-1 to be a cluster blue strag- 
gler, and the star is much too bright to be a horizontal branch or 
blue straggler star in the SMC. An explanation for this object 
remains elusive. 

All of the other blue variables are uncrowded in the Chandra 
image and consistent with non-detection, and their HST time se- 
ries show no evidence for flickering. As with 1V36 this argues 
against them being CVs. Here, we examine possible alternative 
explanations for these blue variables: (a) background variable 
stars from the SMC (e.g. RR Lyraes), (b) detached WD bina- 
ries, (c) very low accretion rate CVs, Uke those discussed in 
Townsley & Bildsten (2002), and (d) exotic collision products. 
A combination of these may also apply. 

Considering possibiUty (a), we plot two SMC RR Lyraes (la- 
beled as Rl and R2) in the WF3 CMDs of Paper I. These have 
periods of 15.2 hr and 8.7 hr and were independently discov- 
ered by two of the authors (Edmonds and Gilliland) using a 
ground-based variability study of a 47 Tuc field off the core 
(14.5' X 14.5'; for background see SiUs et al. 2000). The RR 
Lyrae light curves have a distinctive asymmetrical appearance 
that is clearly different from that of the blue variables. Also, 
only 3V07 hes reasonably close to the expected position of the 
horizontal branch instability strip. A second possibility is that 
the variables are pulsating blue stragglers in the SMC. These 
generally have shorter periods and smaller amplitudes than RR 
Lyraes, for example the range of periods for the main pulsa- 
tion modes of the 6 known SX Phe stars in 47 Tuc are 0.7-2.4 
hr and the range of amplitudes is 0.006-0.085 mags (Gilliland 
et al. 1998). So, while the amplitudes overlap, most of the 
blue variables have periods that are too long. Also, only 3V07 
appears to lie reasonably close to the expected blue straggler 
sequence. Consistent with these negative conclusions, we note 
that the number of blue SMC stars above the MSTO in our field 
is likely to be very small judging by the detection of only 4 
such stars in the study of a WFPC2 47 Tuc field by Zoccali et 
al. (2001). Proper motion studies with separate F300W epochs 



8 



Edmonds et al. 



may help confirm the blue variables as 47 Tuc members. 

PossibiHty (b) is attractive because Fig. 10b shows that sev- 
eral of the blue variables have periods that are reasonably close 
to, but longer than, the values given by equation 1 for Roche 
lobe filling. However, only 2V30 has colors that are obviously 
consistent with those of a MS star/carbon-oxygen (CO) WD bi- 
nary, and the long period for this system of 1.18 or 2.36 days is 
difficult to understand as either ellipsoidal variations or heating 
effects. The colors of 2V08 are very similar to those of a bright 
WD, so any MS companion is likely to be extremely faint. The 
other objects with all three colors available (1V36, 3V07, 4V05 
and the possible ID for W71), have colors, compared to the MS, 
that are bluer inV-I than inU-V, behavior that is difficult to 
reconcile with CO WD-MS star binaries. However, note from 
the color-color CMD given in Paper 1 that 3 V07 lies very close 
to the He WD companion to the MSP 47 Tuc U, suggesting that 
3V07 may also be a He WD. This does not obviously have an 
MSP companion, because of the lack of detected X-ray emis- 
sion, but it could have a WD companion. 

Case (c) appears to be an unlikely explanation because the 
stars are too bright to be consistent with the Townsley & BUd- 
sten (2002) models, and their colors are not obviously consis- 
tent with CO WDs. Without having a satisfactory explanation 
for most of these stars, we turn to possibihty (d) that at least 
some of them represent an exotic collision product formed near 
the center of 47 Tuc. Mathieu et al. (2002) point out, in their 
discussion of red stragglers in M67, that stellar encounters are 
common in star clusters and that it would not be surprising to 
discover products of such encounters (especially binaries) that 
run counter to standard evolutionary theory of single stars. Such 
expectations apply even more so to 47 Tuc. 

4. DISCUSSION 



4.1. Cataclysmic variables 

We have discovered (or confirmed) optical counterparts for 
22 Chandra sources that are CV candidates as summarized in 
Tables 3 and 4 of Paper 1. We have included V3 in this list, 
but, as noted above, it may be a qLMXB. Of these 22 CV can- 
didates, definite variability of some type has been seen for all 
of them except W33opt, W45opt, W49opt, W70opt, W82opt and 
W85opt. The excellent astrometric match between W45 and 
W45opt, the tiny chance of this being a chance coincidence 
(0.012%) and the relatively bright, hard nature of the X-ray 
source makes this ID secure. The other candidates have larger 
possibilities of being chance coincidences (0.8%, 1.0%, 1.4%, 
2.3% and 0.2%) and so without having independent information 
such as variability or an Ha excess there remains the possibil- 
ity that one or even more of these 5 IDs is not real. Also, W70, 
W82 and W85 are relatively soft soures (Fig. 22) and there- 
fore one or two of them may be MSPs (see § 4.3). Among the 
marginal candidates we believe that W140opt and W55opt have 
the greatest chances of being CVs, as explained in Paper I. 

The photometric properties of these stars were sunnmarized 
in Paper I. The other properties of these systems are summa- 
rized as: (1) the Fx /Fop, values are higher, on average, than 
those of all classes of CV in the field sample of VBR97, partly 
because of higher than average X-ray luminosities and partly 
because of faint optical magnitudes; (2) the periods, where 
available, are generally consistent with expectations for Roche 
lobe filling objects; (3) the radial distributions are consistent 

From § 3.3, it is unlikely that qLMXBs, with their lia$e Fx / Fopt values, are si] 
data, only V3 is a reasonable qLMXB candidate 



with those of the cluster MSPs and the brightest MS-MS bina- 
ries, and therefore with WD-MS star binaries; (4) three of the 
CVs are echpsing (W8, W15 and AKO 9), several show signifi- 
cant flickering (VI, V2, V3, & W25opt), two of them have previ- 
ously been seen in outburst (V2 and AKO 9) and several others 
show large amphtude variations (e.g. W51opt, V3, W56opt) that 
may also be signs of outbursts. 

The period/My data presented in § 2.2.3 and the faint optical 
magnitudes presented in § 3.2, plus the apparent lack of non- 
magnetic novalike systems suggest that the 47 Tuc CVs are 
dominated by low accretion rate CVs, i.e. DNe'*. This result 
is not necessarily inconsistent with the lack of outbursts seen 
for 47 Tuc CVs, since some DNe (such as SU UMa systems) 
can have quite long recurrence times (Warner 1995). Shara et 
al. (1996) mention this possibility for 47 Tuc. However, the 
outbursts seen for V2 and AKO 9, and the possible outbursts 
shown in Paper I (the GO-7503 variability) do not appear to be 
as dramatic as those seen in SU UMa systems, and the dis- 
tributions of the 47 Tuc CVs and field SU UMas appear to be 
incompatible. Also, all of the 47 Tuc CVs except for W21opt are 
found above the period gap, while nearly all known SU UMa 
systems are found below the period gap. Therefore, we believe 
it is unlikely that a significant fraction of the 47 Tuc CVs are 
SU UMa systems. 

Selection effects in the field will almost certainly conspire 
against the detection of low accretion rate CVs (Warner 1995). 
For example, nova-like systems are relaitvely bright and blue 
in the optical. Conversely, SU UMa systems like WZ Sge 
are extremely faint at optical and X-ray wavelengths, and have 
very long recurrence times for outbursts. Also, CVs found in 
the field have a range of non-uniform selection criteria. Al- 
though CVs with very low accretion rates (like those modeled in 
Townsley & Bildsten 2002) will be difficult to detect in 47 Tuc, 
the sample of cluster CVs is likely to be relatively complete 
down to X-ray luminosities of 3-5xlO^''erg s~^ Therefore, 
given (1) the depth of our observations, (2) the uniform search 
methods used, and (3) that the 47 Tuc CVs are effectively all at 
the same distance, it may be possible that our 47 Tuc CVs are 
a more representative sample of CVs than the field objects of 
VBR97, and that they therefore show an expected bias towards 
lower accretion rate systems. 

Is the lower metalhcity of 47 Tuc compared to those found 
in field (Population I) CVs likely to result in lower accretion 
rates? Stehle, Kolb, & Ritter (1997) have studied the long-term 
evolution of CVs with low metallicity secondaries (Z=10~^) and 
have shown that such systems have a smaller period gap and a 
sUghtly higher mass transfer rate than CVs where the secondary 
has a solar composition. This trend therefore does not explain 
the low apparent accretion rates for the 47 Tuc CVs. Population 
synthesis modeling by Stehle, Kolb, & Ritter (1997) shows that 
most Population II CVs should be old enough to have already 
evolved beyond the period minimum, and will be extremely 
faint. However, this age effect does not explain the low accre- 
tion rates for the 47 Tuc CVs, since the periods for the cluster 
CVs, where known, are nearly all above the period gap (and 
most of the CVs with unknown periods should also be found 
above the period gap). 

A potential problem with the low accretion rate interpreta- 
tion is that no known class of field CV has both and My 
distributions that are similar to those of the 47 Tuc CVs. As 

gnificantly biasing our results, since based on X-ray spectral information aadFx/Fopt 
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noted above, the X-ray luminosities of the 47 Tuc CVs are 
much higher than those of the U Gem systems, and are con- 
sistent with those of DQ Her systems. So, while the faint op- 
tical magnitudes imply that they have relatively low accretion 
rates, the high X-ray luminosities might imply that they have 
relatively high accretion rates. Only one of these possibilities, 
at most, can be correct. Since My is clearly a better indica- 
tor of accretion rate than Lj, for non-magnetic CVs (VBR97), 
and My is also known to strongly depend on accretion rate for 
magnetic systems (e.g. Hessman, Gansicke, & Mattel 2000 and 
Yi 1994), we speculate that Lx for the 47 Tuc CVs is a poor 
guide to their accretion rates whether they are magnetic or non- 
magnetic systems. We note that two systems with presumably 
relatively low accretion rates are V2 (a DN) and W56opt, a pos- 
sible DN. These are two of the brightest CVs in X-rays, with 
Lx= 6.3 X 10^'erg s~' and 1.5 x lO^^erg s"^ respectively, hinting 
that high X-ray luminosities are compatible with low accretion 
rates. We encourage theoretical work on globular cluster CVs 
to try to explain these results, and caution that the low accretion 
rate interpretation given here is not secure until the combined 
X-ray and optical luminosities are understood. 

Given the relatively high X-ray luminosities of the brightest 
CVs in 47 Tuc and the observation that about 25% of field CVs 
are magnetic (Wickramasinghe & Ferrario 2000), is it possible 
that most of the brighter CVs are DQ Her systems? Solid evi- 
dence exists for DQ Her-like behavior in VI (GHEOla), and the 
strong resemblance between the spectra of AKO 9 and GK Per 
suggests that AKO 9 is also a DQ Her system (Knigge et al. 
2001). Unfortunately, these two systems represent just a small 
fraction of the total CV population in 47 Tuc, and the period vs 
My and Fx/Fopt data presented in § 2.2.3 and § 3.2 shows that 
the 47 Tuc CVs are statistically different from DQ Her systems 
found in the field (field DQ Hers have, on average, significantly 
brighter optical magnitudes; Fig. 21b). This contrasts with the 
similarity between the distributions of the 47 Tuc CVs and 
the field DQ Her systems (see Fig. 21a). 

The and My values for the 47 Tuc CVs are consistent 
with the data presented for NGC 6397 by Cool et al. (1998), 
Edmonds et al. (1999) and GHEOla. The values of the 
NGC 6397 CVs are high (like those in 47 Tuc), since 4 out 
of 9 systems have Lj. > lO^' ** erg s"'(GHE01b). Optically the 
NGC 6397 CVs are faint, despite being expected to have even 
brighter secondaries than 47 Tuc (because of the low metalUc- 
ity, [Fe/H]=-1.95; Harris 1996). The absolute magnitudes of 
the 4 relatively bright CVs from Cool et al. (1998) and Ed- 
monds et al. (1999) range between 5.95 and 8.81. 

Further information about the NGC 6397 CVs are available 
because good quality optical spectra have been obtained. Ar- 
guments have been made that the CVs in NGC 6397 may be 
dominated by DQ Hers (Grindlay et al. 1995 and Edmonds et 
al. 1999), based on the observation that 3 of the 4 NGC 6397 
CVs with measured optical spectra have relatively strong He II 
4686A lines, like those found in field DQ Hers. Figure 8 of 
Edmonds et al. (1999) shows an analysis of the absolute mag- 
nitudes of the accretion disks [Mi^(disk)] for DQ Her systems 
and the ratio between the continuum levels of H/3 and Ha 
[C(H/?)/C(Ha)] for different classes of CV. For field DQ Her 
systems, linear relationships are found between My (disk) and 
the He II 4686A to H/3 equivalent width ratio (Hell/H/3), and 
a similar relationship is found between [C(H/3)/C(Ha)] and 
HeII/H/3 (brighter and bluer disks, equivalent to higher ac- 
cretion rates, correlate with larger Hell/H/J values). The 4 



NGC 6397 CVs discussed in Edmonds et al. (1999) have 
My (disk), [C(H/3)/C(Ha)] and HeII/H/3 values that are con- 
sistent with the field DQ Her systems, but intriguingly, the 3 
NGC 6397 CVs with moderate HeH/H^ values are all found in 
the low accretion rate part of the figure, with low Mv(disk) and 
low C(H/3)/C(Ha) values. Therefore, the low accretion rates 
apparently found for the 47 Tuc CVs, combined with the hints 
of magnetic behavior in a couple of systems, may be consis- 
tent with the NGC 6397 results. Optical spectra of the fainter 
CVs in NGC 6397 and the large population of CVs in 47 Tuc 
are needed to test whether they also have relatively strong He II 
4686A (useful spectra will be difficult to obtain because the 
CVs are mostly faint and crowded). Searches for rapid WD spin 
in the optical and X-rays is also needed. Perhaps these obser- 
vations will identify a new class of low accretion rate, magnetic 
CVs in globular clusters. 

One final possibility worth considering is that somehow the 
different formation mechanism for globular cluster CVs com- 
pared to field CVs is responsible for differences in their respec- 
tive behavior. To test this hypothesis we note that CVs in open 
clusters are much more likely to be formed from primordial 
binaries than the 47 Tuc CVs. Interestingly, the two CVs dis- 
covered in the open cluster NGC 6791 by Kaluzny et al. (1997) 
are both much bluer (in quiescence) in the V vs V -/ CMD than 
most of the 47 Tuc and NGC 6397 CVs (see Edmonds et al. 
1999), suggesting brighter accretion disks and higher accretion 
rates than the 47 Tuc and NGC 6397 systems. Indeed, one of 
the NGC 6791 CVs appears to be either a UX UMa or Z Cam 
system, with a relatively high accretion rate. Only one other CV 
is known in an open cluster, a faint AM Her system (EU Cancri) 
in M67 (GilHland et al. 1991, see Fig. 10; Belloni, Verbunt, & 
Mathieu 1998). Clearly, any two randomly selected CVs from 
47 Tuc would not look Uke the NGC 6791 CVs, but obviously 
the statistics are very poor. Larger samples of CVs in rich open 
clusters and low density globular clusters are needed to see if 
this effect is significant, or just a statistical anomaly. 

4.2. Active binaries: comparisonwith field studies 

We have a total of 29 likely active binaries in 47 Tuc. Twenty 
seven of these show statistically significant, mostly periodic, 
variability, and most of them are found on, or slightly above the 
MS or subgiant ridge-line except for a handful of red stragglers 
or red straggler candidates. The total sample of active binaries 
in 47 Tuc above our X-ray detection threshold will inevitably 
be larger, since a sensitive variabiUty study was only possible 
with the GO-8267 data. 

A key question is whether the active binaries presented here 
have similar properties to those of BY Dras and RS C Vns found 
in the field. Here, we present a brief comparison with the 
ROSAT All-Sky Survey results of Dempsey et al. (1997). Fig- 
ure 4 of Dempsey et al. (1997) shows that the X-ray lumi- 
nosities of RS CVns and BY Dras are quite different from each 
other, with median X-ray luminosities of ~ 1 .6 x 10^' erg s"' for 
the BY Dras and ~ 4 x lO^^erg s"^ for the RS CVns and max- 
imum luminosities of 2.5 x lO-'^erg s"' for the BY Dras and 
~ 1 X lO^^erg s~' for the RS CVns. These X-ray luminosities 
are the quiescent values, as Dempsey et al. (1997) observed 
flares on all of their targets during their survey and removed 
data points with enhanced count levels. 

The red stragglers (W3opt, W43opt and W72opt), the red strag- 
gler candidates (W4opt, W14opt and W38opt) and the stars near 
the MSTO like W9opt and W75opt generally have X-ray lu- 
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minosities (see GHEOla) that are in good agreement with the 
Dempsey et al. (1997) data if the primary stars are subgiants 
and these systems are RS CVns. The two red stragglers (and 
possible sub-subgiants) in M67 both have L^= 7.3 x lO^^erg s"' 
(Belloni, Verbunt, & Mathieu 1998), in good agreement with 
the X-ray luminosities of W3opt and W72opt. Also, the rela- 
tively bright V mags (and in some cases obviously red colors) 
of the X-ray bright BY Dras (Fig. 16d) are consistent with the 
hypothesis that a reasonable number of these objects have at 
least partially evolved primaries and hence are RS CVns. 

Among the active binaries with optically fainter IDs the 
source W47 is probably returning from a flare, and X-ray vari- 
ability in other sources such as W18, W41 and W94 may also 
be signs of flares, explaining the relatively high luminosities 
even though these systems are likely BY Dras. The luminosi- 
ties for the opticaUy faint and apparently non-X-ray variable 
systems like W26, W59 and W66 are obviously more extreme. 
Some of them could have X-ray variability that is intrinsically 
quite large, but is not detected because of the faintness of the 
sources. Another possibility is that they are simply outliers in a 
large population of objects having a distribution that is similar 
to those found in the field. AGBOl point out that large numbers 
of BY Dras may have been missed in their study due to in- 
completeness, since the binary frequency determined from the 
BY Dra sample (with periods that are mostly a few days long) 
is a factor of 17 smaller than that calculated from the shorter 
period eclipsing and contact binaries. Assuming that the binary 
frequency is the same among the BY Dras in the AGBOl study 
as it is among the short period ecUpsing and contact binaries, 
and scaling from the sample of 55 AGBOl binaries with V > 19 
this would imply a population of about a thousand BY Dras 
in 47 Tuc. Since Figure 4 of Dempsey et al. (1997) shows 
that ~ 15% of the field BY Dras have > 10^°erg s"', we 
would therefore have a deficit of X-ray detected BY Dras with 
Lx > lO^^erg s~^ , rather than a surplus. Possible reasons for this 
deficit are that a reasonable fraction of the binaries with periods 
longer than about 1-2 days have been destroyed by interactions, 
or that there is an anti-correlation between Lj^ and period (see 
below). 

The plot of X-ray color vs V in Fig. 22 shows that most of the 
optically faint active binaries are relatively soft sources, consis- 
tent with the expected soft X-ray spectra for these objects in 
quiescence (Dempsey et al. 1993 and Dempsey et al. 1997). 
The source W47 is a relatively hard source and is variable, and 
flares are expected to harden the X-ray spectrum (as noted in 
Paper 1 W64 may show long-term variability). 

If the bright active binaries have a substantial fraction of 
RS CVns and the faint systems are all BY Dras then there may 
be an anti-correlation between Lj and V (ie the optically fainter 
active binaries being associated with fainter X-ray sources). 
We found no statistically significant evidence for such an anti- 
correlation, either with or without the X-ray variables (by con- 
trast the CVs show weak anti-correlations between and V, 
especially when AKO 9 is removed from the sample). With 
a much larger number of (mostly undetected) BY Dras than 
RS CVns it is much more likely that a few very luminous 
BY Dras are detected rather than a few very luminous RS CVns, 
possibly explaining the lack of correlation noted above. 

In agreement with the trends shown in Fig. 17d, there is a sig- 
nificant anti-correlation between period and V mag for the ac- 
tive binaries. The linear correlations between period and V are 
-0.38 and -0.48 when using two different tests (the 'kendll' 



and 'pearsn' subroutines) from Numerical Recipes (Press et al. 
1992), with chance probabilities of only 0.4% and 1%. Presum- 
ably the main cause of this anti-correlation is that the brighter 
active binaries contain a number of systems with evolved stars 
and relatively long periods. However, why have only the short 
period BY Dras (the faint active binaries) been detected in X- 
rays? According to Dempsey et al. (1997), there is a fairly 
weak correlation between X-ray flux (F^) as measured at the 
Earth and rotation period (Pro,) of: F^^ for BY Dras 

(a stronger correlation exists for RS CVns). Assuming that 
Lx= Pi^t'^^ for the active binaries, then the median period of 
0.4 days for the faint active binaries (the BY Dras) and 1.56 
days for the faint AGBOl binaries means that L, would be only 
17% lower for the AGBOl binaries. This is probably not a large 
enough difference to explain the very different period distribu- 
tions, although a steeper slope (but still remaining within the 
1-fT error limit) could easily produce a 30-40% difference be- 
tween the X-ray luminosity of the faint active binaries and the 
faint AGBOl binaries. We also note that the period distribution 
of the Dempsey et al. (1997) sample is very different from ours, 
since only 3/29 of the Dempsey et al. (1997) BY Dras with 
known periods have periods < 1 day, but 15/27 of the 47 Tuc 
active binaries with known periods have periods < 1 day (our 
photometrically selected sample is biased towards short peri- 
ods, explaining part of this effect). With this different sample 
of periods (with a much stronger bias towards short periods) it 
is not surprising that we find a different Fi/period relationship. 

In an important sense our sample of active binaries in 47 Tuc 
(at fixed metalUcity and age) is significantly cleaner than the 
Dempsey et al. (1997) sample containing an inhomogeneous 
sample of stars with a range of ages and metallicities. There- 
fore, the examination of trends such as the dependence of F^ on 
period may be optimally performed with cluster samples (we 
defer this study to a future pubHcation foUowing analysis of 
deeper Chandra observations obtained in late 2002). 

4.3. Are some of the active binaries MSPs? 

Here, we examine the possibiUty that some of the active bi- 
nary candidates could be MSPs. Such a population would be 
very interesting for dynamical reasons (see § 1), and it has a 
direct impact on the number of MSPs in 47 Tuc estimated from 
the Chandra data. The large range in the estimated number 
of 47 Tuc MSPs (35-90) given by Grindlay et al. (2002) is 
dominated by assumptions about the nature of the active binary 
population (see also § 4.4). 

We begin by noting that the eclipsing binary X-ray sources 
(W12opt, W92opt, W137opt, & W182opt) are unlikely to be MSPs 
(or CVs) because two eclipses would not be expected. Also, 
the X-ray sources in the AGBOl sample classified as W UMas 
(W41opt, W47opt, & W163opt) have periods and colors that obey 
the Rucinski relationship for such systems (AGBOl), therefore 
making the contact binary explanation a more natural one, al- 
though here we do not rule out an MSP explanation. 

The radial distributions of the active binaries presented in 
§3.1 offer powerful constraints on this issue. There are striking 
similarities between the radial distributions of the active bina- 
ries and the AGBOl binaries (with the active binaries removed) 
in both the bright and faint samples (see § 3.1 and Fig. 17). The 
similarity in the bright samples suggests that the bright active 
binaries are not dominated by MSPs with 'normal' MS com- 
panions, because then mass segregation would mean that these 
~2.15-2.25M0 objects (1.4M0 neutron star + 0.75-0.85Mq 
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MS star) should be much more concentrated towards the center 
of the cluster than the lower mass objects (the radio detected 
MSPs and the CVs). There is a chance that some of the bright 
active binaries could be MSPs with low mass (< O.2M0) MS 
star companions that have been heated by X-rays so that they 
appear optically like MSTO stars or red stragglers, as may be 
occurring for 6397-A (Ferraro et al. 2001). In principle the 
radial distribution of these objects could then mimic those of 
MSTO binaries, however, this would require a coincidence and 
would not explain an apparent absence of MSP binaries with 
MS star companions having masses close to the turn-off value 
of ~O.85M0. A second possibility is that binary-binary and 
binary-single star interactions resulting in dynamical kicks may 
cause the average MSP-MS star binary to be further away from 
the center of the cluster than expected for mass segregation of 
such objects. However, this would require another coinciden- 
tal similarity between radial distributions. It would also be in- 
consistent with the difference noted in § 3.1 between the radial 
distributions of the bright and the faint active binaries. This dif- 
ference probably would not have been seen if MSPs dominate 
both the faint and the bright groups because then the percent- 
age difference in mass between the two groups would have been 
much smaller (as shown in Fig. 17 the faint active binaries are 
unUkely to contain a significant number of MSPs based on their 
radial distribution). 

When these radial distribution results are combined with the 
evidence that (1) the faint active binaries fall above the MS 
ridge-line just like the AGBOl sample of binaries (see Paper 1), 
and (2) that the X-ray luminosities are consistent with RS CVn 
and BY Dra systems (§ 4.2), we conclude that the active bina- 
ries are indeed dominated by MS -MS binaries. 

Is this finding consistent with the properties of the known 
MSP population? Of the 16 MSPs with timing positions (Freire 
et al. 2001) only half of them are binaries and of these only 5 
have companion masses likely to be above the minimum value 
(O.O85M0) for being on the MS (the other 3 have masses of 
~ 0.02-0.03Mq). It is likely that these 5 MSPs all have He 
WD companions based on the radio data (Camilo et al. 2000), 
e.g. 47 Tuc U has an //5r-identified He WD companion (Ed- 
monds et al. 2001) while a faint, blue star and possible ID has 
been found for 47 Tuc T and faint limits can be set in the optical 
for 47 Tuc H. Only one MS companion was found (47 Tuc W; 
Edmonds et al. 2002b) in the full sample of 20 MSPs from 
Camilo et al. (2000). Therefore, if the sample of 16 MSPs 
is representative of the possible much larger sample of MSPs 
(Camilo et al. 2000) then optical companions are rare and MS 
companions even more so. If a large population of MSPs with 
MS companions does exist, and we are detecting many of them 
in X-rays, then the MSPs would have to be intrinsically radio- 
faint to avoid any of them being detected by Parkes. However, 
the relatively flat correlation between X-ray and radio luminos- 
ity shown by Grindlay et al. (2002) impUes that the faintest 
MSPs in the radio are less faint as X-ray sources. 

Given the small number of detections of cluster MSPs with 
likely non-degenerate companions (two; 47 Tuc W and 6397-A) 
it is not yet possible to conclude whether such MSPs are gener- 
ally relatively bright in X-rays but faint in the radio. The MSP 
6397-A would have been easily detected by Parkes if it were 
in 47 Tuc, given the radio luminosity of 5 mJy kpc^ quoted by 
D'Amico et al. (2001). The radio luminosities of the 47 Tuc 
MSPs range between 0.81 and 10.9 mJy kpc^, with 12 out of 
14 of them having radio luminosities less than that of 6397-A. 



The MSP 47 Tuc W, however, is one of the faintest detected 
MSPs, and W34, if it is an MSP, is below the detection limit of 
Camilo et al. (2000). The X-ray luminosities for 6397-A and 
47 Tuc W are slightly higher than the most luminous 'normal' 
MSP in 47 Tuc. 

Without having a deeper sample of MSPs, and detailed radial 
velocity information for the active binaries it is very difficult to 
rule out the possibility that any individual, non-eclipsing, opti- 
cal variable has an MSP companion, especially in view of the 
cases of W29opt and 6397-A (and possibly also W34opt). The 
former object is very crowded optically and was only discov- 
ered using difference image analysis (Edmonds et al. 2002b), 
but was found on close study to have blue colors. Therefore, 
W29opt would have been classified as a CV instead of an active 
binary if not for the remarkable period and phase coincidence 
with 47 Tuc W discovered by Edmonds et al. (2002b), perhaps 
implying that one or two of the faint CV candidates could be 
an MSP. None of the active binary candidates show blue colors, 
and they also show no hint of the large amplitude variability 
seen in either W29opt or W34opt, despite in some cases having 
periods that are not much longer than those of W29opt (suggest- 
ing similar levels of variability caused by irradiation from an 
MSP, if present). 

One potential source of MSPs masquerading as active bina- 
ries are the red stragglers. For example, the companion to the 
MSP 6397-A was previously identified as a BY Dra by Taylor 
et al. (2001) on the basis of Ha emission, a CMD position ly- 
ing above the MS and optical variability (this star is clearly a 
red straggler). However, it was subsequently shown by Ferraro 
et al. (2001) to have an MSP companion based on an accu- 
rate radio timing position for the MSP. Naturally, it is possi- 
ble for a binary system to be both an active binary and to have 
an MSP companion, since secondaries in short period binaries 
Uke 6397-A are undoubtedly tidally locked and rapidly rotat- 
ing, and hence may have substantial chromospheric activity (as 
suggested for 6397-A by Orosz & van Kerkwijk 2002). 

Although the 47 Tuc red stragglers and 6397-A have similar 
optical properties, there are some key differences between their 
X-ray properties. While 6397-A is a reasonably hard source (as 
is 47 Tuc W), the two red stragglers W3 and W72 both have 
very soft X-ray colors and the two red straggler candidates W4 
and W14 also have soft colors (see Fig. 22 and Fig. 24). Of the 
red stragglers, only W43 (in GO-7503) has a hard color. Gen- 
erally, if most of the active binary candidates are MSPs and are 
like 47 Tuc W and 6397-A, then they should mostly be hard 
sources. Yet, of the GO-8267 active binaries with X-ray color 
determinations from Grindlay et al. (2002) only 4 (not includ- 
ing W47opt) are hard while 9 (not including W4) are soft. 

The leading candidates for MSPs hidden within our active 
binary sample include W43, a faint, moderately hard, variable 
X-ray source (like 6397-A) with an optical ID that falls in a 
similar region of the optical CMD to 6397-A, as noted above. 
No clear optical variability is present, with our limited quality 
time series. We note that while W43 is clearly harder than the 
other red stragglers, this could be because of flaring activity, 
explaining the variability. Both W23 and W64 are moderately 
hard sources that were not observed to vary during the Chandra 
observation of GHEOla, and are brighter than all of the other 
active binaries (except W47). However, as explained earlier, 
W64 may show long-term variabihty. 

The leading candidates for identifications of other 47 Tuc U- 
like MSPs (not yet detected in the radio), are the weak, soft, 
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non-variable X-ray sources, with faint, blue, optical compan- 
ions that do not show long-term variability. These are W82 and 
W85 (see Fig. 22 and Fig. 24), with plausible optical IDs (lying 
outside the GO-8267 FoV) without obvious long-term variabil- 
ity. They are beyond the faint limit of the GO-7503 data for 
short-term variability studies. The ID for W70 shows no evi- 
dence for periodic variation in our time series and only marginal 
evidence for non-periodic variations, but we do not yet have 
long-term variability information for this object. These sources 
could be CVs, since AM Her systems for example are known 
to be soft X-ray sources (VBR97). Our high quality Ha obser- 
vations with HST/ ACS may help distinguish between these two 
possibilities, since He WDs should show broad Ha absorption 
lines, and CVs strong Ha emission lines, as observed for the 
CVs in NGC 6397 (Grindlay et al. 1995 & Cool et al. 1998) 
and NGC 6752 (Bailyn et al. 1996). 

4.4. Revised estimate of total number ofMSPs 

Finally, we examine the soft, faint, apparently non-variable 
X-ray sources that have no plausible (optically variable or pho- 
tometrically unusual) counterparts (see Fig. 24b). Grindlay et 
al. (2002) listed the sources in this category, noting in partic- 
ular the subset of these sources that fall in the GO-8267 FoV: 
W4, W5, W6, W24, W28, & W98 (we now add the source W71 
to this hst). As noted by Grindlay et al. (2002) these sources 
include the leading candidates to be MSPs. By correcting for 
the incomplete spatial coverage of the WFPC2 GO-8267 FoV, 
and for incompleteness in the detection of radio-detected MSPs 
as soft X-ray sources, Grindlay et al. (2002) estimated the pres- 
ence of at least 19 MSPs with > 10^°erg s"' that do not have 
radio counterparts with timing positions, in addition to the 16 
MSPs with such positions. This estimate obviously increases 
if some of the active binaries identified with soft sources are 
instead MSPs, but as argued above, few are expected. Also, 
the estimate by Grindlay et al. (2002) imphcitly assumes that 
2/9 of the MSPs have X-ray colors < 1.0, so a small fraction 
of misidentified 'active binaries' with harder X-ray colors (like 
W9 or W59) that are really MSPs, do not change these statis- 
tics. 

Based on a study of the optical data, is it reasonable that the 
7 soft sources just discussed are all MSPs, rather than CVs or 
active binaries in some cases? First, as already noted, W4 is a 
candidate red straggler and W24 may be an RS CVn, and both 
of them have higher X-ray count levels than the brightest de- 
tected MSP (except for the hard X-ray source 47 Tuc W). Also, 
W98 is located close to several giant stars so it may also be 
an RS CVn. To estimate how many active binaries, or perhaps 
CVs, may have missed detection because of crowding, we have 
examined the stellar count levels present in the GO-8267 data 
around the positions of the 7 soft sources. Figure 23 shows 
the mean stellar counts in an annulus (between 0."9 and 1 "35) 
around each: (1) CV, (2) bright active binary, (3) faint active 
binary (4) MSP, (5) unidentified X-ray source, and (6) candi- 
date MSP (the 7 sources discussed above). Optical positions 
are used for (1), (2), (3) and (6), radio positions for (4) and 
X-ray positions for (5). The median count level for the can- 
didate MSPs is higher than any of these other groups, and is 
2.0 times larger than the median count level for the MSPs (the 
two distributions differ at the 98.8% level using the KS-test). 
This suggests that some of the 7 soft sources are active bina- 
ries (or CVs) that have been missed because of crowding. We 
tentatively classify W4, W24 and W98 as active binaries, based 



on their promixity to bright stars, and classify the 4 remaining 
sources as MSPs. This imphes a population of 13 extra MSPs 
for a total of 29 MSPs with > lO^" erg s'^. 

4.5. Summary of optical identifications 

Incorporating the optical identifications reported here and 
in Paper I (plus radio/MSP identifications from GHEOla and 
Grindlay et al. 2002) and just considering the GO-8267 FoV 
alone (where our high quality time series coverage is spatially 
complete) a significant fraction of the 78 sources are active 
binaries (34.6%) compared to 19.2% for CVs and 9.0% for 
MSPs, with the source W46/X7 being a qLMXB (1.3%). These 
sources, where X-ray colors are available, are labeled in the X- 
ray CMD shown in Fig. 24. The remaining 27 (34.6%) sources 
in the GO-8267 FoV have uncertain identifications, or have no 
plausible ID. Of these sources, we use the X-ray luminosity and 
color, or the plausible match-up with bright, blue or marginally 
variable stars to tentatively identify WIO, W16, W17, W20, 
W32, W35, W37 and W140 as CVs, and W4, W24, W37, 
W71, W93, W141 and W168 as active binaries leaving 13 faint 
sources that we identify as possible MSPs (see Table 3). This 
leaves us with final fractions of 42.3% active binaries, 29.5% 
CVs and 26.9% MSPs (plus the single qLMXB), although some 
of the 13 weak X-ray sources could be CVs or active binaries, 
lowering the MSP number The CV fraction is very similar 
to that calculated by GHEOla, but the active binary fraction is 
significantly higher than the 15% estimated by GHEOla. Cor- 
respondingly, the MSP estimate, and the one given in § 4.4, are 
smaller than the value of ~50 MSPs originally presented for the 
108 sources in GHEOla (that estimate was based on incomplete 
analysis of the optical data). 

4.6. Conclusion and Prospects 

Based on the period/My analysis and the My distribution pre- 
sented here, the 47 Tuc CVs appear to have lower accretion 
rates than field CVs systems found above the period gap. The- 
oretical work on CV evolution predicts that the fainter, shorter 
period CVs expected to exist in 47 Tuc (but not observed be- 
cause of crowding) should have even lower accretion rates than 
the bright ones that we do observe (see Di Stefano & Rappaport 
1994 and Townsley & Bildsten 2002). The suggestive evidence 
for low accretion rates and DQ Her type behavior in CVs in 
NGC 6397, combined with: (1) the evidence for magnetic be- 
havior in a couple of the 47 Tuc CVs (VI and AKO 9) and (2) 
the high (DQ Her-like) X-ray luminosities of the 47 Tuc CVs 
implies that magnetic activity may play an important role in 
globular cluster CVs. However, more work is needed to under- 
stand the unusual My and distributions for the CVs in both 
47 Tuc and NGC 6397. 

Perhaps the most unexpected result, when compared with 
other clusters, is that the active binaries outnumber the CVs. 
In NGC 6397, the detected CVs outnumber the active binaries 
by a factor of ~2-3 (based on deeper Chandra data than ob- 
tained for 47 Tuc), but significant depletion of the MS binary 
population in NGC 6397 is likely to have occurred during and 
after core collapse. In NGC 6752, the detected CVs outnum- 
ber the active binaries by a factor of ~4-10, but the NGC 6752 
Chandra data is not as deep as the 47 Tuc data, and several 
unidentified soft sources may still be identified with active bi- 
naries, given optical data of higher quality. 

With this large population of detected active binaries in 
47 Tuc, many of them soft sources, we find support for the 
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lower range given by Grindlay et al. (2002) for the number 
of MSPs in 47 Tuc with > lO^^erg s"' (~30^0). It is possi- 
ble that a few of the active binaries may really be MSPs, and we 
continue to search for new objects that are similar to 47 Tuc U 
or 47 Tuc W, but, given the spatial, photometric. X-ray spec- 
tral and timing information available for the observed active 
binaries, a large population of MSPs masquerading as active 
binaries is unlikely. The lower estimated number of MSPs in 
47 Tuc means that the ratio between the number of MSPs and 
the number of CVs is smaller, and closer to the corresponding 
value found in NGC 6397 by GHEOlb. 

A comparison of Figures 18 and 19a shows that we are sen- 
sitive to detection of most classes of CV except for the faintest 
DNe (SU UMas), and possibly also double degenerate systems 
(where the statistics are poor). There are good prospects for 
finding at least some of the fainter CV population not yet de- 
tected here by using the much deeper X-ray limits set by the 
300 ks ACIS-S observation. Although the gains will be lim- 
ited near cluster center because of X-ray crowding, the faintest 
CVs should have low-mass companions of ~ 0.1 -O.3M0, and 
with WD masses of ~O.55M0 the total system masses should 



be less than or equal to the MSTO mass. Therefore, many of 
these objects should be found in relatively uncrowded regions 
more than 20-30" away from the cluster center. Unfortunately, 
optical identifications will be difficult or impossible for many 
of these objects because of crowding and the high background 
levels present in the HST images, a penalty for studying objects 
in this massive, concentrated globular cluster. 

The 300 ks foUowup Chandra observations will provide 
much better quality spectral, timing and positional information, 
especially for the faintest sources discussed here, and these data 
should be very useful for distinguishing between MSPs and ac- 
tive binaries. Many new, faint sources will also likely be dis- 
covered, dominated by MSPs and active binaries. New radio 
observations are obviously needed to search for the ~200 un- 
detected MSPs believed to exist in 47 Tuc (Camilo et al. 2000). 

We acknowledge comments on the CVs from Brian Warner, 
and a number of extremely helpful comments from the ref- 
eree. Christian Knigge. This work was supported in part by 
STScI grants GO-8267.01-97A (PDE and RLG) and HST-AR- 
09199.01-A(PDE). 
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Fig. 1. — Phase plots and Lomb-Scargle power spectra in the V and / band for Wl (a CV). Individual points and phase bins (with 3-it errors) are shown for 
the phase plots and a sinusoid fit is overplotted. The orbital period (and half of it) are labeled on the power spectra, and the HST orbital frequency is shown. The 
horizontal dotted line shows the power level corresponding to a false-alarm probability of 1 x 10"^. 
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Fig. 2. — Phase plots and power spectra in the V and / band for the eclipsing CVs W8opt and W15opt. Note, for both systems, the deeper eclipses in the V band 
compared to the / band because the blue accretion disk is being eclipsed. These CVs both have very hard X-ray spectra, likely caused by significant X-ray absorption 
by the accretion disk. 
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Fig. 5. — Power spectra in the V and /band for VI, V2, V3 and W45opt- Possible periods are labeled (for elUpsoidal variations the periods will be twice the values 
shown). The X-ray counterpart to V3 (W27) has a period of 3.83 hr (GHEOla). 
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Fig. 9. — Finding charts (in F300W) for the CVs or CV candidates found in the GO-7503 FoV that show detectable or significant long-term variability. Data from 
GO-5912 (epoch 1), GO-6467 (epoch 2) and GO-7503 (epoch 3) are shown (the appropriate epoch number is shown in parentheses). The solid lines show the 4~a 
error circles (except in the cases of W56 and W120 where 8-cr error circles are shown for clarity) and the dashed circles encircle the IDs (with radius O'.'l for PCI 
and Of' 1 for the WF chips). 
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Fig. 1 0. — A plot of orbital period vs My for 47 Tuc CVs and CV candidates (Fig. lOa) and blue variables (Fig. 10b). In Fig. 10a two CVs from each of the 
globular clusters NGC 6397 (CVl and CV6; Grindlay et al. 2001b and Kaluzny & Thompson 2002) and NGC 6752 (Bailyn et al. 1996) are shown plus a CV from 
M5 (Neill et al. 2002) and M67 (Gilliland et al. 1991). Small filled triangles are used to plot the non 47 Tuc CVs. The periods for the blue variables have been 
assumed to be twice those found from the pealc in the power spectrum. The straight line shows equation 1 from Warner (1995). 



22 



Edmonds et al. 




Optical Counterparts to 47 Tuc X-ray sources. II 



23 




Fig. 12. — Cumulative plots of the difference between My predicted by equation 1 (using the binary period) and the measured My for 47 Tuc CVs (labeled as 
'Tuc CVs) compared to field CVs from VBR97. This difference is an estimate of the brightness of the accretion disk and hot spot, plus (for faint systems) the WD, 
where higher accretion luminosities are towards the right of the figure. The secondaries in the 47 Tuc systems will be brighter than solar metallicity secondaries, 
and so the 47 Tuc accretion luminosities wiU be even fainter than shown by these figures. The 'DNe' distribution is for all DN systems. Fig. 12a shows all systems 
with periods less than 10 hr and Fig. 12b shows systems with periods less than 10 hr but above the CV period gap (this shows systems where the contribution from 
the white dwarf is likely to be smaller). 
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Fig. 14. — Time series for a sample of active binaries with long periods. These are a mixture of red stragglers (W3opt and W72opt) and red straggler candidates 
(W14opt), binaries and likely BY Dra variables located near or above the MSTO (W9opt, W73opt and W75opt) and a variable located well above the MS ridge-line 
(W69opt). Short horizontal segments in the time series show various gaps in the time series. 
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Fig. 16. — Cumulative radial distributions for several classes of 47 Tuc binary: CVs, MSPs, binaries detected by AGBOl ('binaries'), the subset of these binaries 
('xABs') detected in X-rays and the blue variables (Fig. 16a); The V distribution of this same subset of binaries, where the MSP list includes only 47 Tuc U, 
47 Tuc W (W29), W34 and the possible ID for 47 Tuc T (Fig. 16b); The period distribution for these binaries, where the periods are derived from AGBOl, Camilo et 
al. (2000) and this work (Fig. 16c); The V distribution for the AGBOl binaries classified as BY Dra systems and the subset of these systems that are X-ray sources 
(Fig. 16d). 
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Fig. 17. — Radial distributions for bright (V < 19; Fig 17a) and faint (V > 19; Fig 17b) subsets of the AGBOl binaries and the active binaries. Fig. 17c shows the 
radial distributions of the bright and faint active binaries, the total stellar population (mostly single stars) in the quoted V ranges, and the unidentified X-ray sources 
in the GO-8267 field of view; Fig. 17d shows the period distributions for the bright and faint AGBOl binaries and the active binaries. 
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Fig. 19. — Plots of Fx/ Fopi vsMy for the 47 Tuc optical identifications of Chandra sources. The symbols used as the same as Fig. 8, except that arrows are shown 
for limits in Fx (blue variables) or My (faint CVs in GO-7503 and 47 Tuc T). 
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Fig. 20. — Cumulative plots of Fx /Fopt for the 47 Tuc CVs compared to the 47 Tuc active binaries and MSPs (Fig. 20a) and various classes of field CV from the 
work of VBR97 (Fig. 20b). For the non-magnetic field CVs the accretion rates decrease from left to right. These results imply that the non-magnetic CVs in 47 Tuc 
have low accretion rates, although the Fx /Fopt values are biased by the unusually high Lx values for the 47 Tuc CVs. 




Fig. 21. — Cumulative plots of (Fig. 21a) and My (Fig. 21b) for the 47 Tuc CVs compared to the field CVs from VBR97. Upper limits have been included 
in the field Lx distributions to increase the sample size. Note the similarity between the distribution of the 47 Tuc CVs and the field DQ Her systems, and the 
similarity between the My distribution of the 47 Tuc CVs and the field U Gem and AM Her systems. 
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Fig. 22. — X-ray color plotted against V . For the X-ray color we use the definition of GHEOla and the selection criteria of Grindlay et al. (2002) for faint sources. 
Fig. 22a shows the CVs, MSPs, the qLMXB X5 and CV candidates, and Fig. 22b shows the active binaries. 
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mean annulus counts 

Fig. 23. — Cumulative distribution of mean count values (stellar and background) for annular regions in the GO-8267 V-band image centered on CVs, bright 
(y < 19) active binaries, faint (V > 19) active binaries, MSPs, unidentified X-ray sources and candidate MSPs (soft X-ray sources without optical IDs). The optical 
positions are used for the CVs and active binaries, the optically-corrected radio positions are used for the MSPs (except for 47 Tuc W where the optical position is 
used), and the optically-corrected X-ray positions are used for the unidentified sources and the candidate MSPs. 
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Fig. 24. — X-ray CMD plotting Lx in the 0.5-2.5 keV band versus X-ray color (as defined by GHEOla). Fig. 24a shows the optical identifications of CVs, MSPs 
and active binaries in the GO-8267 dataset, the radio/X-ray identifications of MSPs without optical IDs (47 Tuc D, E, G/I, H, N; Camilo et al. 2000; Freire et al. 
2001) and X7 (W46) an X-ray identified qLMXB (GHEOla and Heinke et al. 2003) without an optical ID (Edmonds et al. 2002a). Fig. 24b shows sources with 
X-ray colors but lacking optical or radio identifications (some of these have marginal optical IDs). 



Optical Counterparts to 47 Tuc X-ray sources. II 



Table 1 

Fx/Fopt VALUES FOR 47 TUC SOURCES AND FIELD CVS 



Source type 


Number 


(median) 


xABs 


29 


0.0059 


MSPs 


4 


0.46 


CVs 


17 


1.18 


Field CVs" 


UX Umas 


9 


0.0072 


Z Cams 


9 


0.043 


AM Hers 


9 


0.22 


DQ Hers 


8 


0.23 


U Gems 


13 


0.39 


SUUMas 


18 


0.73 



"CVs from the study of Verbunt et al. 
(1997) 



Table 2 

ks-test probabilities for fx / f,,,,, distributions being different 



Source SUUMas UGems DQHers AM Hers ZCams UXUMas CVs MSPs 



U Gems" 


77.39 












DQ Hers" 


86.13 


13.40 










AM Hers" 


97.06 


42.68 


0.82 








ZCams" 


99.84 


99.13 


87.52 


92.22 






UXUmas" 


100.00 


99.99 


99.26 


99.65 


98.13 




CVs" 


97.72 


98.20 


93.36 


98.44 


99.93 


100.00 


MSPs" 


41.28 


1.13 


1.43 


4.57 


90.46 


99.11 


xABs" 


100.00 


100.00 


99.73 


99.93 


99.74 


8.17 



"CVs from Verbunt et al. (1997) 
"Optically identified sources in 47 Tuc 



Edmonds et al. 



Table 3 

Chandra SOURCES with marginal or no //5r counterparts in the GO-8267 data 

Source Chip* Crowding level Plausible 
W# ID" 



zu 






i_ V ; 


24 




B 


AB? 


28 




A 


MSP? 


31 




B 


MSP? 


32 




D 


CV? 


35 




C 


CV? 


37 




B 


CV? 


39 




B 


47 Tuc O" 


46 




B 


qLMXB 
MSP?' 


96 




A 


97 
98 




C 
D 


MSP?' 
AB? 


141 




C 


AB?' 


4 


2 


D 


AB? 


5 


2 


D 


MSP? 


10 


2 


D 


CV? 


13 


2 


C 


47 Tuc N'' 


16 


2 


D 


CV? 


17 


2 


D 


CV? 


19 


2 


C 


47 Tuc 0,1" 


67 


2 


C 


47 Tuc D'' 


74 


2 


C 


47 Tuc H'' 


91 


2 


A 


MSP?' 


95 


2 


B 


MSP?' 


142 


2 


D 


MSP?' 


168 


2 


D 


AB?' 


115 


3 


C 


MSP?' 


6 


4 


D 


MSP?' 


7 


4 


D 


47 Tuc E" 


71 


4 


C 


MSP? 


93 


4 


C 


AB?' 


99 


4 


A 


MSP?' 


140 


4 


A 


CV?' 



Note. — Grade for crowding level: (A) Source in 
uncrowded region; (B) Source in moderately crowded 
region; (C) Source in extremely crowded region; (D) 
Source in region affected by saturation 

''Plausible source identification based on X-ray spec- 
tral information and nearby objects in the optical 

•"MSP from Camilo et al. (2000) 

'^Source does not formally have an X-ray color, using 
the criterion of Grindlay et al. (2002) 



